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A theoretical calculation is presented to describe the confined motion of transmembrane molecules in cell 
membranes. The study is analytic, based on Master equations for the probability of the molecules moving as 
random walkers, and leads to explicit usable solutions including expressions for the molecular mean square 
displacement and effective diffusion constants. One outcome is a detailed understanding of the dependence of 
the time variation of the mean square displacement on the initial placement of the molecule within the con- 
fined region. How to use the calculations is illustrated by extracting (confinement) compartment sizes from 
experimentally reported published observations from single particle tracking experiments on the diffusion of 
gold-tagged G-protein coupled yu-opioid receptors in the normal rat kidney cell membrane, and by further com- 
paring the analytical results to observations on the diffusion of phospholipids, also in normal rat kidney cells. 
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I. INTRODUCTION 

Much activity has recently centered around the biophysics 
of cell membranes. This interest stems from the variety of 
processes in which the membrane plays a key role, among 
them cell shaping and movement cell division |2|, sig- 
nal transduction |3|, and molecule trafficking |4J. This latter 
process is of fundamental value for the regulation of the local- 
ization, assembly, and aggregation of molecules within, and in 
the vicinity of, the plasma membrane, functionalities that are 
all linked to a healthy cell existence Q. Understanding the 
motion of membrane-associated molecules is, thus, of direct 
interest and great relevance. 

Recent observations of lateral movements of molecules on 
the surface of the cell lEllgllTlElilllOllIl | have led to the sug- 
gestion that the moving (transmembrane) molecules are con- 
fined within certain regions of the cell membrane. This con- 
finement is ascribed [ 1 1 1 to collisions of membrane molecules 
protruding into the cytoplasm, the substance that mainly fills 
the cell interior, with the cytoskeleton |12|, a filament-like 
structure present in the interior of the cell. The collisions are 
expected to reduce the movements of the molecules and ef- 
fectively confine their motion on the surface of the cell, as if 
the plasma membrane were compartmentalized. Within each 
compartment, the molecules are envisaged as moving freely, 
their motion being hampered as they traverse adjacent com- 
partments. As the actin filament that forms the compartment 
boundary dissociates due to thermal fluctuations, a nearby 
transmolecule with sufficient kinetic energy may overcome 
the barrier potential and hop to the adjacent compartment. 
These arguments are the basis of the so-called "membrane- 
skeleton fence model" for temporary "coralling" of the trans- 
membrane molecules lfT3l . Further slowing eff'ects are also 
possible due to "transmembrane-protein pickets" anchored on 
the membrane skeleton fence [11 J, in agreement with the lat- 
est research findings ||T41 . suggesting that the lipid bilayer 
plasma membrane is more mosaic than previously thought 
i\5\. A recent publication by Morone et al. |fT6l on imag- 
ing the cell surface by using electron tomography techniques 
has provided further support for these models. 



Observations IfTTl [T3l [TTl of molecular movements in the 
plasma membrane of various mammalian cell types point to 
compartment dimensions between 30 and 240 nm and an aver- 
age hop rate between compartments ranging from 1 to 17 ms. 
Whereas diffusion constants of these molecules are known to 
be about 10 yum^i"' in the absence of compartments, the above 
observations lead to the extraction of effective diffusion con- 
stants 5 to 50 times smaller for the motion associated with 
hops between compartments 1, 11 J . These experimental results 
have been analyzed IfTTl in the past in terms of the model of 
Powles et al. [IS], wherein the molecule is looked upon as 
a random walker moving in an infinite 1-D space with peri- 
odically arranged semipermeable barriers. The short-coming 
of the theory delineated in ref. [T8 1 is that the space-time de- 
pendence of the probability distribution of the random walker 
is obtained in terms of an unwieldy infinite series of terms 
which is difficult to handle. Furthermore the mean square dis- 
placement, which is the quantity of direct comparison to the 
experiment [ll j, is not derivable analytically from ref. ITSl . 

Against this backdrop, which consists of experiments 
exhibiting non-standard (meaning non-free) diffusion of 
molecules on a cell membrane, and a theory which is not easy 
to validate, use, or manipulate, because of its involved ex- 
pressions, we present a new analysis for hindered molecular 
diffusion on cell membranes. It is based on calculations car- 
ried out by one of the present authors, but not published, on 
the formally (but not physically) similar problem of Frenkel 
exciton transport f2A \ in deuterated molecular crystals. We 
will see that our theory results in easily usable expressions for 
the molecular mean square displacement and for the effective 
diffusion constant. We will show explicitly how to use our ex- 
pressions by deducing reasonable values for the compartment 
size from experiments on ceU membranes. 

The paper is organized as follows: The model and its for- 
mal solution are in section [H] the continuum limit leading to 
the effective difiiision constant is in section III and explicit 



usable expressions for the time dependence of the molecular 
mean square displacement are in section IV A comparison 



to the experiments of Kusumi and collaborators iTTJI is given 
in section |V] and concluding remarks constitute section |VI] A 
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few calculational details may be found in the Appendix. 



II. THE MODEL AND ITS FORMAL SOLUTION IN 
DISCRETE SPACE 

In constructing our model we follow Powles et al. fTSl 
in considering a random walker in an infinite 1-D space but 
we study a discrete chain of sites with lattice constant a. Al- 
though, for the sake of simplicity, we display the results only 
in 1-D as in earlier work 1 18 1, generalization to higher dimen- 
sions (e.g., 2-D) is straightforward at least in principle. The 
molecule, whose probability of occupation of the m-th site of 
the chain at time t is Pmit), hops with nearest neighbor rates 
F within a compartment of H +1 average number of sites and 
with a lower rate / from the end site of the compartment to 
the first site of the adjacent compartment. For simplicity, H is 
taken to be even here, with the site at the center of one of 
the compartments. The equation of motion is therefore 



dPM 
dt 



if m is within a compartment 
dPrif) 



dt 



f [Pr^iit) - P,{t)] + F [P,_i(f) - P,{t)] .. 



for the end site of one compartment, and 

dPr^lit) 



dt 



: / [Pr{t) - Pr^m + F [Pr+2{t) " ^r+l(f)] : 



for the first site of the next compartment. Here r and r+ 1 mark 
the site location of, respectively, the left and right ends of each 
compartment boundary. The index r equals H/2+(H+l)l with 
I an integer running from -oo to +oo. 

With A-F-f, the three equations can be combined: 



dP,n{t) 

dt 



= F[P,„+i(f)+Pm-l(f)-2P„,(f)] 



-A J] [Pr+lit) - P,.(t)] - (5„,,,+ i) (1) 

r 

where the primed summation is over the left ends r of all bar- 
riers and the 5's are Kronecker deltas. Equation Q is formally 
identical to an equation studied in the past (see ref. ^\ and 
references therein) to describe exciton processes in a molecu- 
lar crystal. Using tildes to denote Laplace transforms and e to 
denote the Laplace variable, it is possible to solve Eq. ([TJ as 
(see Appendix ) 

X\^^^r(e)-^,n-r-x{e)\, (2) 

where (f) - ^,n-n (0 is the probability propagator for the 
system without barriers (A - 0), i.e. the probability that 
the molecule is at site m at time t if it was at site n at time 



in the barrier-less system. The dependence of the prop- 
agator on merely the difference in the site locations stems 
from translational invariance. By t]„ (f) is meant the solution 
Yin "^m-n (0 Pn (0) of the homogcncous (barrier-less) problem 
for the given initial probability distribution P„ (0). 

Information about the location of the barriers is in the cru- 
cial function yu(e), which, (see Appendix), for periodic place- 
ment of barriers, is given by 



tanh (^/2) 



tanh(^(//-H l)/2) 



1 



(3) 



where cosh^ =1-1- el2F. This closed form for Jiie) is a 
consequence of the fact that the propagators for the barrier- 
less chain with nearest-neighbor rates, i.e., the system obey- 
ing Eq. ([T]) with A = 0, are products of modified Bessel 
functions and exponentials in the time domain, specifically 
*?/(?) - Ii{2Ft)e^^^' , the consequent Laplace transforms be- 



ing 



-fl'l 



2F sinh ^ 



Since the ?;'s are easily obtainable (being solutions of the 
homogeneous equation, as explained above), in light of Eq. 
Q, Eq. (|2]i constitutes the full solution for the probabilities. 
A quantity of direct contact to experiments on the cell mem- 
brane is the mean square displacement. Because of the pres- 
ence of barriers, one has to take care in defining this quantity, 
focusing in particular on its dependence of where the molecule 
starts in relation to the barriers. Translational invariance of the 
compartments allows one to focus on any one compartment 
without loss of generality. We consider thus, as the most gen- 
eral case for a fully site-localized initial condition, the particle 
starting at the p''^ site inside the central compartment so that 
P„(0) - 6„ p, p being any integer in the interval [-H/2, H/2]. 
In this case, rj,„ - ^P,,,-;,. Then the mean square displace- 
ment for a locaUzed initial condition, defined as the average 
of (m - pf-. 



• pfPr, 



(t). 



and labeled by the starting site p, is given in the Laplace do- 
main by: 

yp(e)^Yj^m-pf^,„^p(e) 

m 

f 

-A J] pM) Yj.m - pf ^m-r{e) - ^,„-r-x{e)\ , (4) 

r m 

and has obvious dependence on the starting site p. Here pr, 
defined and evaluated in the Appendix, is simply the differ- 
ence - P,-. Since = (;7t_j r)^ + 2r(m - r) + P', we 
can write, using the fact that Ym = 1 /f and that symmetry 
considerations dictate that the first moment of the propagators 
vanishes. 



Yjm - pf \Y^-Ae) - 'Vn,-,-x{e)\ = - 



2r-2p + \ 
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obtaining the mean square displacement as the sum of its 
counterpart for the barrier-less case and a correction which 
is proportional to the difference of the rates A-F-f: 



(5) 



(6) 



By defining a function gp{t) with the Laplace transform 

and writing (pp (t) - 5{t)- gp(t), we have 

yp(t) = 2F f dt' f if>p{t")dt" 
Jo Jo 

where (pp{t) may be called the memory function [2T]. 

One of the advantages of expressing the mean square dis- 
placement in terms of the memory is that, for any starting site 
p, the effective hopping rate F^ff at large times can be easily 

r+co 

calculated as F^ff - F dt(pp(t) and is given by 



eff 



F<Pp(0)^f 



H+ 1 



\+{flF)H 



(7) 



For illustrative purposes we will show a few examples of 
gp{t) here, labeling them explicitly as gH,p(f) to show the com- 
partment size as well. When H -2, i.e., when a compartment 
is composed of only three sites, P{e) - 1 + e/2F, Q(e) = 1 
and 



g2,o(t) = 



-3Fi 



g2,±i(t) = A| 



3F e 



-3Fr 



2 f-F 



(F^-Ff/2 + f)e 
F(f-F) 



-{F+2f)t 



2F I 



(9) 



When H - A each compartment possesses 5 sites, P{e) - 
1 + 2elF + e^/2F^, Q(e) = 2 + e/F and we obtain 



g4,o(0 = lOF 



3^ sinh ( VSFf /2) 



sinh(V5-4/+4(/)V^/2) 



(10) 



That ^(0), and consequently (pp(0), is independentt of p will 
be clear below. Equation (j7|l shows that if / = F, there is 
no confinement effect and F^ff - F. If the intercompartmen- 
tal motion is absent, i.e., if / = 0, the molecule cannot es- 
cape and Feff - 0. For this case the mean square dispacement 
saturates to within the compartment. For small intercompart- 
mental escape, i.e., if / << F, specifically if we can neglect 
( f/F)H « 1, then F,ff = f(H+l). This means that the 
molecule jumps large distances of compartment size at the 
lower rate /. We also see, the prediction of little interest to the 
cell membrane problem but interesting in other contexts, that 
if on arrival at a barrier the molecule is whisked away into the 
next compartment, specifically if / >> F/H, there is a small 
enhancement of the motion rate: F^ff = F (l + ^)which can- 
not exceed a factor of 3/2. 

In the Appendix it is shown that the evaluation of Eq. (|6]l 



gives 



A 



1 



FF{e) + (el2F + 2flF)Q{e) 
P{e) + {el2F)Q{e) 
^ ' P(e) + (2 + e/2F)Q(e) 

2p sinh(^p) tanh(^/2) 



cosh(^p) 



(8) 



where P(e) = cosh(^///2) and Qie) = sinh(^///2)/ sinh^. 
Note that 'gp is symmetric in the sign of the initial site p as 
expected and that for e = 0, and therefore for ^ = 0, it be- 
comes independent of p as stated above. 



Returning to the original 'gp in Eq. ([8]), and averaging it over 
all initial p within the compartment with equal weight, gives 
a /:> -independent barrier contribution to the memory: 



Hf2 



1 



F(H+l) 



p=-H/2 

P(6) + (e/2F + 2)Q(6) 
P(€) + (e/2F + 2f/F) Q(e) 



(11) 



Note that g(t) — > as either f — > oo or f — > 0, since eg{e) 
as e — > and e — > oo. 



III. CONTINUUM LIMIT 



It is straightforward to compare our analytical result with 
corresponding Monte-Carlo simulations of random walks 
with barriers. We performed such simulations by following 
standard procedures f25l in a 1-D lattice with N{H + 1) sites 
with large enough so that the molecule (random walker) 
never reaches the boundaries during one run of the simulation. 
In Fig. [T] the analytical results and the results of the simula- 
tion are shown for H - A. We have averaged over 20000 
different trajectories for each case to create the final results. 
Time is displayed in units of 



4 




written as 



5 10 

time (in units of 1/F) 
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FIG. 1: Illustrative comparison of our analytic result for the mean 
square displacement with Monte-Carlo simulations, showing excel- 
lent agreement. Parameters are / = 0.01, F = 2 H = 4. Solid 
curves show the result of the simulations, averaged over 20000 runs. 
Dashed lines are the analytical results. They are almost indistinguish- 
able from one another. The inset shows the behavior at short times. 



{(x - xof) (/) - 2Dt - 2£' j| dt' j| g',^ (f") dt". (14) 

Notice that^'^(e) maintains the limiting properties of ^(e), i.e., 
e^'(e) — > as e — > and e — > oo indicating that (f) — > as 
f — > oo and f — > 0. As in the discrete case, at large times it is 
possible to calculate the diffusion coefficient, 

Jr»oo 
fif/[(5(f)-/(f)]=-D[l-?^(0)], 




which turns out to be 



D.ff = 



D 



1 + 



(15) 



This result, which appears as the continuum limit of our more 
general discrete expressions, is equivalent to that obtained in 
ref. IfTSl . 



While the above clear coincidence of the simulation results 
with our analytic solutions establishes confidence in the anal- 
ysis, the comparison has been made only for small H for sim- 
plicity. Practical application of our analytic result in the cell 
membrane context requires an infinite number of sites within a 
compartment, i.e., the continuum limit, because the molecules 
diffuse, rather than hop, from location to location. Therefore, 
we now proceed to the continuum, assuming that, as the lat- 
tice constant a — > 0, the products Ha, pa, and fa tend to 
the respective constants L, xq, and S!, but that F tends to in- 
finity an order of magnitude faster so that Fa- tends to (the 
diffusion constant) D. Surely, xq is the initial position of the 
particle inside the central compartment and 2! ID corresponds 
to the permeability of a partially permeable barrier, indicated 
in ref. ifTSl as V. Similar continuum limits may be found in 
ref. ||26l. By retaining only the nonvanishing terms in Eq. ([8|, 
the continuum Umit is given by 



Us)- 



s cosh s + sinh s 



coth s cosh ( s 
L 



2x0 . , /2xo 

sinh s 

L \ L 



(12) 



where the dimensionless quantity s - ^ yf^ is proportional 
to the square root of the Laplace variable, the subscript jcq de- 
notes the initial location on the continuum, and the superscript 
c is a reminder of the fact that we are treating the continuum 
limit. 



The continuum limit of Eq. (Ill becomes (this corresponds 
to an average over xq) 



sinh s 



s cosh s + 



sinh s 



(13) 



In terms of the inverse Laplace transform g'^ (t) of Eq. ( 12 1, 
the mean square displacement in the continuum limit can be 



IV. EXPLICIT USABLE EXPRESSIONS FOR THE MEAN 
SQUARE DISPLACEMENT 

For most of the rest of the paper we will examine only the 
continuum limit expressions because they allow direct com- 
parison to experiment in the specific problem of molecular 
motion in cell membranes. Other systems such as those in- 
volving the hopping of excitations in a molecular solid I2TII 
are better treated with the help of the more general discrete 
case expressions we have obtained above. We will now use the 
notation y^^ to mean {{x - xq)^) and describe how to take the 



inverse Laplace transform of Eq.( 12 1 to get the mean square 



displacement in the time domain, with the help of Eq. ( 14 1. 
Note that the mean square displacement in the Laplace domain 
can be written as 



Deff 



1 



1 



8D \ cosh s 

1 tanh s 



i{s;xo,L), (16) 



where 



^{s; xq, L) - coth s cosh | 



2x0 . . /2xo \ 

sinh s\ 

L \ L 



(17) 



Average of Eq. ( 16 1 over all initial locations xq gives 

-L/2 



y(s) 



- r 



L 

8D 
1 



dx(){(x - xoYKs) 



L/2 
1 



tanh s 
tanh^ s 



tanh.+ ^i' 



(18) 
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We notice that in Eqs. (16 1 and (fTSll, the mean square dis 



placement expression in the Laplace domain appears decom- 
posed naturally as a sum of two parts. The first term repre- 
sents the mean square displacement for the case of impenetra- 
ble barriers. It can be inverted easily to give a simple analytic 
expression in the time domain: 



I? 
2 



32 y 



11=0 



1)4 



(19) 



This result for the mean square displacement of a totally con- 
fined particle and related expressions for the autocorrelation 
function of the displacement agree with published expressions 
that have appeared earlier in the study of nuclear magnetic 
resonance microscopy and of animal motion in home ranges 

MM- 

Many problems in physics are tackled in terms of an ex- 
actly soluble system to which a perturbation is added. The 
choice of the partition of the system into an unperturbed and a 
perturbed part is obviously never unique and is itself an inter- 
esting aspect of theoretical study. This happens, for instance, 
in quasiparticle transport in molecular crystals and aggregates 
ETl wherein a localized and a delocalized part of the system 
represent two such possible choices. We see here that our 
present system displays such a choice in that one may regard 
the unperturbed part of our random walk problem as being ei- 
ther the walk on the free chain perturbed by the presence of 
the multiple barriers, or as the confined motion of the walker 
within a single well with two impenetrable barriers (/ - 0) 
perturbed by a leakage to the rest of the chain. In the first 
case, represented by Eq. (|5]l, the unperturbed system is rep- 
resented by free motion on a chain. In the second, with the 
partition described by Eq. ( [T8] l, it is Eq.( 19 1 that describes the 
unperturbed system. 

We now calculate the inverse Laplace transform of the full 
mean square displacement, focusing directly on^(e)/e^ since 
it is the double time integral of g(t) that is required in the cal- 
culation. Here, and henceforth in this section, we drop the 
superscript c on ^ to avoid clutter and do not display the ex- 
plicit dependence on the starting location xq. The well-known 
identity |24J 



X-i{/(aVi)}(f) 



1 



f 

Jo 



duf{u)ue 



(20) 



has found use earlier in problems similar to the present study 
f27J . It allows one to invert Laplace transforms of functions 
of the square root of the Laplace variable if the corresponding 
transforms of just the variable are known. 
In the time domain, we have 



y(f) = 2Df-2D r dt' r dt" 
Jo Jo 

deg{e)e"\ 



(21) 



where C is the Bromwich contour enclosing all the singular- 
ities of 'gie) (note that all of the singular points of Eq. ( 12 1 



are on the imaginary axis). After performing some simple (al- 
though tedious) algebraic manipulations, employing the nota- 
tion y-Dgff I { D - Dgff ), a-2xo/L, T-ADtlL?, and using 
Eq. ([20|, we get 



y{r) = y (T - icri (r) - y(T2{T) - acT,{T))) , (22) 



where 



1 C cosh(Q'i)e^ 

o-i(t) = — (p as — r— — , 

2m Jc smh s 



o"2(t) 



0-3(7) 



1 



ds- 



cosh(Q's)e'"^ 



2m Jq {s cosh s + y sinh s) 



2m ^ 



sinh(Q;s)e ' 



ds- 

c i (i cosh i H- y sinh i) 



(23) 



Therefore cri, cr^, 0-3 will be equal to the sum of the residues 



of the integrands in Eqs. (23 1. Let 



wi(s,t) = 



cosh(Q's)e'' ^ 

sinh s 
cosh(a^)e''''" 



W3(i, r) - 



s^{s cosh s + y sinh s) ' 
sinh(a5)e'''" 



s(cosh s + y sinh s). 



(24) 



Then, 



- y ^ Res[w2, i„) - ff ^ Res[wi, s„]^, (25) 



where s„\ are the solutions of -s - ytanh,?. At s - 0, 
wi and W2 have poles of order 3 and W3 has a pole of or- 
der 2. The corresponding residues are t + a^/2 - 1/6, 
(6T(y + 1) + 3a^(y + 1) - (y + 3))/(6(y + if) and a/iy + 1) 
for wi, W2 and W3 respectively. Note that each of the functions 
wi, W2 and W3 has infinitely many poles besides the pole at 
s=0. The rest of the poles of wi are all simple poles (as can 
easily be seen by Taylor-expanding sinh s) located at s-imn 
where m - +1, ±2, ±3... . For W2 and W3, we cannot find the 
exact locations of the poles analytically because of the need to 
solve a transcendental equation. Although we do not have an 
analytical expression for the location of these poles, we know 
that they are all simple poles because the first derivative of 
s cosh s + y sinh s does not vanish at s„. Therefore, excluding 
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the pole at s„ = in the n-summations below, 



cri(T) = T + - 1/6 + y lim 



cosh(Q'i)e*''" 
cosh s 



(26) 



0"2(t) 



6T(r + 1) + 3a;% + 1) - (r + 3) 
6(r+ 1)2 

Zcosh(Q'i)e'''" 
lim , 
^ s^s„ sinh s + (1 + y) cosh s) 

a sinh(ci's)e*''" 

<T3(t) = — + > lim — — — — 

(y + 1) ^ s^s„ sis sinh 5 + (1 + 



• (1 + 'y)coshi) 



(27) 



(28) 



The roots of l/sinh(s) can be found analytically, and cri(T) 
can be expressed as 



^i(t)= £ (-l)"" 



J cos(amn)e 



(29) 



the right hand side being simply related to the r-integral of 
elliptic theta functions. 

Because we cannot solve -s - ytanhi exactly, it is not 
possible to write analytical expressions for cr2{T) and cr3(r). 
However, we can find the roots s„ of that equation numerically 



with high precision, and evaluate the sums in Eqs. (27i and 



(28 I to get cr2(T) and cr3(r). In doing this, we use the bisection 
method to find the first few thousands of the roots of cr2{T) 
with an accuracy of 10 Because of the factor e""^ in Eqs. 



(27 1 and (28 1, the sums converge very fast as i„'s are purely 
imaginary. Note that lim,^ cosh(a5) and limj^ ,,^ sinh(Q'.s) lie 
in [-1, 1] and the magnitude of lim^^j^ ((1+y) cosh sinh .?) 
tends to oo with increasing lm(5„). As a result, it is possible to 
obtain accurate results without summing over a large number 
of residues. A similar procedure can be followed f29^| to find 



the inverse Laplace transform of Eq. ( 18 i 



We have displayed several of the details of the inversion 
procedure we have used because of its involved nature and 
because it is not often that (even partial) inversions resulting 
in explicitly usable time domain expressions are possible in 
similar problems. In Fig. |2]we show the dependence, on the 
(dimensionless) time t, of the normalized instantaneous diffu- 
sion coefficient. 



Pit) 
D 



1 d 
2Ddt 



y(t) 



2 d 
l?d^ 



y(r). 



(30) 



V. COMPARISON WITH EXPERIMENT 

In comparing our theoretical results with experimental ob- 
servations in the cell membrane field, we used two sets of 
data by digitizing two time plots of the molecular mean 
square displacement available in the published literature. The 
first is Fig 4b (left) in ref. |19| which is about the diffu- 
sion of a gold-tagged G-protein coupled yU-opioid receptor; 




FIG. 2: Time dependence of the mean square displacement and the 
instantaneous diffusion coefficient for various magnitudes of the con- 
finement effect. Plotted as a function of time is the normalized y(t) 
(left) and D(t)/D (right), averaged over all initial locations in the 
compartment, for different values of Dgff/D: 0.001 (solid line), 0.01 
(dashed line), 0. 1 (dash-dotted line), 0.5 (dotted line). All quantities 
considered are made dimensionless appropriately. 



the second is Fig 2b in ref. which is about the diffu- 
sion of a phospholipid molecule (l,2-dioleoylsn-glycero-3- 
phosphoethanolamine). Both sets of observations refer to nor- 
mal rat kidney cells. The observations were made at a time 
resolution of ISfis which is high enough to capture the change 
in the diffusion constant at times long and short compared to 
expected values of L-/4D. 



We first fitted our y vs f curve, given by Eq. (21 1, to the 
data in Fig 4b of ref. |19|, both to the x and y components 
of the mean square displacement, and extracted the compart- 
ment sizes through standard procedures. It is not possible to 
make accurate deductions because the data in the literature are 
available only for single trajectories rather than for ensemble 
averages for a large number of them. Within this limitation, 
we have been able to conclude by the application of our theory 
to the data that the linear extent of the confining compartment 
should lie between 250mot and AlQnm. Our conclusion is in 
agreement with the distribution of compartment sizes given in 
Fig 4d of ref. lfT9ll and generally with values discussed in ref. 

ma. 



A different system we examined is that involving the diffu- 
sion of smaller transmembrane molecules than proteins: phos- 
pholipids 1 20 1. By evaluating Eq. (21 
deduced values for D, De/f and L, 



using experimentally 



we verified that conclu- 
sions derived in ref. [201 are in agreement with theory. This is 
displayed in Fig. |3] To show the agreement clearly, we use for 
the compartment size, L, a substantially different value from 
that deduced in ref. 1201, 230 nm, and point out that the y vs 
t curves found by using Eq. ( [2T| do not then agree with the 
experiment at all. 

As the data corresponds to a single trajectory, it is worth- 
while to explore the effects of the initial position of the 
molecule on the dynamics. In Fig. |4] we show the effect of the 
initial position of the particle on the yvst curve. The param- 
eters D, Deff and L are the same as the ones reported in ref. 
EOl . For illustrative purposes we also show the y vs f curve 
averaged over all initial positions as the dashed line in Fig. |4] 
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FIG. 3: Comparison of diffusion observations of phospholipids re- 
ported in ref. |20il to our theoretical predictions based on 3 disparate 
values of the compartment size L used in our Eq. |2TJ. One of the 
values used is that deduced in ref. f20l, 230 nm. There is very good 
agreement of theory (solid line) and observations as shown. Each 
of the other two values differs by an order of magnitude: 2300 nm 
(dashed line) and 23 nm (dotted line). In both these cases theory 
shows poor agreement with experiment. We have used £)=4.6 yura^i-"' 
and Deff=l.2 fim^s^^ . 
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0.01 



FIG. 4: Theoretical predictions for various initial placements of the 
molecule compared to observations reported in ref. |20|. Circles and 
triangles denote observations. Lines represent our theory. Dashed 
line shows the average of the mean square displacement over all ini- 
tial positions. For the solid lines, the initial position is: xo=0, L/8, 
L/4, 3L/8 and 7L/16, going upwards from the lowest curve in the 
plot. Values of D and D^ff that correspond to the observations are as 
in Fig. [3]and the deduced compartment size is L=230 nm. 



VI. CONCLUDING REMARKS 

The investigation reported in this paper has had two sepa- 
rate aims. One is to develop an efficient and useful formalism 
to describe the motion of a transmembrane molecule on a cel- 
lular membrane following ideas inherent in the current mem- 
brane skeleton fence model of Kusumi and collaborators 11 IJ . 



The other aim is to understand, in general, the problem of a 
random walker moving in a region with confinement. Work- 
ing in a 1-D system for simplicity, we gave the exact solution 
of the probabilities in the Laplace domain for a discrete chain, 
derived from them analytic solutions for the mean square dis- 
placement of the molecule, and introduced memory functions 
to understand the evolution physically. We achieved a detailed 
understanding of the motion arising from repeated encounters 
of the random walker with the barriers and dependence on ini- 
tial placement of the walker. These discrete solutions are of 
use not only to the present problem but also to diverse contexts 
such as the motion of excitations in crystals and molecular ag- 
gregates. 

We took the continuum limit of our results to obtain them 
in terms of parameters that can be measured experimentally in 
the context of the cellular membrane and presented a proce- 
dure to get explicit and easily usable expressions in the time 
domain. While a piece of the Laplace inversion has to be done 
either approximately or numerically, other physically relevant 
parts are implemented fully analytically, the entire calculation 
being quite usable from a practical standpoint. Our expres- 
sions are more useful in this regard than those available ear- 
lier, e.g., in ref |18|. Specifically, our results may be typified 
by the prescription for the instantaneous diffusion coefficient 
normalized to the diffusion constant, i.e., D{t)/D given in Eq. 
(30 1. ExpUcitly, identifying the m- summation | 



2 (-l)'"cos(a7rm)e-'^'""' 



with the elliptic theta function of the fourth kind, '&4, as given 
in Eq. (8.180) of ref. |l28l, we can write. 



y lime'' 



j„,s„#0 



D(r) ^ Deff 
D D 

^COSh(Q'i) 

cosh s 



+ 1 - ■&4 (an/l I inr) 

Del, 



D-D,, 



+ as tanh(a'i) 



D 



D-D,f 



+ s tanh s 



(31) 



Here a, defined earlier as 2xo/L, measures the initial location 
relative to the compartment size, and t = ADtjL^ is a dimen- 
sionless time. At infinite time, t — > oo, since the theta function 
tends to 1 and the summation over the s„'s vanishes because of 
the exponential factor (note that that the s's lie on the imagi- 
nary axis), Dij) becomes the effective diffusive constant £>c//. 
At the initial time, on the other hand, the theta function van- 
ishes and D{t) takes the value D for all initial placements of 
the molecule not at the edge, xq Ljl. The mean square dis- 
placement is obtained trivially from the integral of the right 
hand side of Eq. 



(31 



Of the three terms to be computed in 
Eq. ( [3T| ), the first two are given analytically in our theory and 
the last is obtained numerically. 

Our expressions are characterized by three parameters: L, 
the compartment size, D, the initial diffusion constant (as t 
0), and Deff, the effective diffusion constant (as t oo). We 
showed how our theory can be used directly to deduce the 
compartment size L by fitting the analytic expressions to data 
obtained in single particle tracking experiments li 1 fl [191 . 
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---p=0 
p = +4 




Time (in units of 1/F) 

FIG. 5: Effect of initial placement of the random walker at site p 
within a compartment on the time dependence of the effective trans- 
fer rate which, in the continuum limit, would be proportional to the 
time-dependent effective diffusion coefficient. Each compartment 
has 11 sites, i.e., H = 10, and f/F = 0.01. Averaging equally over 
all initially localized placements results, as shown by the dash-dotted 
line, in a monotonically decreasing transfer rate. However, interest- 
ing structures appear for initial placements at the center of, end of, 
and elsewhere in the compartment as shown. See text for explana- 
tion. 



An interesting outcome of our analysis is the description of 
the explicit influence of the initial placement of the random 
walker. The influence is observed in the time dependence of 
the mean square displacement. The dependence is evidently 
a consequence of the relative location of the barriers and the 
site of initial placement. With the help of the analytical ex- 
pressions we have obtained, we plot in Fig. |5] the instan- 
taneous time derivative of the mean square displacement, a 
quantity that is proportional to the instantaneous diffusion co- 
efficient D{t) in the continuum limit, for several cases of the 
initial placement p of the molecule within the compartment. 
We show the discrete case so that the results can be carried 
over to other physical contexts such as excitation transfer as 
well. We observe that the time derivative of the mean square 
displacement does not simply change monotonically from an 
initial value to a lower final value. Instead, interesting struc- 
tures appear. 

The compartment considered in the plot has a total of 1 1 
sites. For central initial placement {p = 0), we see (this is the 
top curve shown dashed) the appearance of a curious dip at 
the bottom of the curve. The dip means that for some time the 
molecule diffuses slower than what it does eventually via the 
final effective diffusion constant. The situation changes quite 
a bit as we vary the site of initial placement, although the final 
effective diffusion constant always remains unaffected. For 
initial placement at the edge of the compartment, denoted in 
Fig. |5]by the dotted curve {p - +5), the molecule begins with 
an effective transfer rate which is not 2F as in other cases but 
a lower value F + f.lt then rises first, reaches a peak, and then 
decreases to the final effective value. For initial placement 



which is neither at the center of the compartment nor at its 
edge, shown in the plot by the solid curve (p - ±4), we see 
that the effective motion occurs with an initial transfer rate 
2F and decreases, both features being shared with the case of 
central placement, but that there is a subsequent increase to a 
peak and decrease to the final effective value, features shared 
with the case of edge placement. 

It is easy to understand all these features. The molecule, if 
placed centrally within the compartment, tends to move ini- 
tially with the transfer rate or diffusion constant characteris- 
tic of the barrier-less system until it meets the barrier. At this 
point it crosses the barrier more slowly and the effective trans- 
fer rate drops. When the molecule has diffused to the next 
compartment it is outside the immediate influence of the bar- 
rier and the effective diffusion is therefore faster. The com- 
bined effect of repeated free diffusion and baiTier-hindered 
diffusion eventually brings the effective diffusion constant to 
its long time value. If, however, the initial placement of the 
molecule is at the compartment edge, it already begins moving 
with an effective diffusion constant lower than the free value 
because of the immediate effect of the barrier. It diffuses faster 
from then on until other barrier encounters including the one 
at the other edge of the initial compartment decrease the rate 
of diffusion. For intermediate initial placements these effects 
happen one after the other as the molecule encounters first the 
barrier on one side and then the one on the other side of the 
compartment. 

Whether the structures described are visible in a specific 
kind of experiment will obviously depend on the freedom one 
has in initial placement of the random walker. We believe in 
light of the data that we have seen in the cell membrane field 
that the average alone matters and that the average, as shown 
in the plot, shows merely a monotonic decrease of the effective 
diffusion constant. The dips have been seen in the theoretical 
work of Powles et al 1 18 1. However, no explanation has been 
given. Indeed, (see Fig. 3 of ref. |18J and the subsequent 
discussion) the authors of that work referred to the dip as an 
"unexpected minimum" and did not connect it to the dynamics 
we have described above. 

Missing from our present analysis are the effects of disorder 
in the placement and height of the barriers, i.e., in the random 
nature of H and / in real systems. Static disorder effects of 
this nature are being analyzed by us in several different ways, 
including the replacement of the periodic expressions for the 
quantity p used in the present paper by random counterparts 
that take into account probability distributions, and the use of 
effective medium theories. These analyses will be reported in 
a separate publication. Of interest to us are also calculations 
of random walks in the presence of dynamic fences that allow 
us to address properly the case where the moving molecules 
are the lipids while the fences are made by the proteins (as has 
been suggested in the literature |11|, and calculations which 
are valid in the presence of intermolecular interactions. 
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APPENDIX A 

1. Evaluation of 

Exploiting the linearity of Eq. ([T]i, it is possible to write its 
solution in the Laplace domain as 

r 

where we denote the difference between the probability to the 
left and right of the barriers Pr+\{t) - PAt) by prit) and its 
Laplace transform by Pr{e). To find an explicit solution from 
the above equation, we first write its particular cases for m - s 
and m - s + \ where s and s + \ mark, respectively, the sites 
to the left and to the right of a barrier: 

Ps^xie) = 77.+i(e) - A ^ p,{e) ^,-,-^i{€) - ^.^-Ae)] , 

r 

f 

PAe) = l]s(e) - A 2 p,{e) ^ .-A^) - ^ .-r-Ae)\ 

r 

(Al) 



and then subtract the second expression in Eq. ( Al i from the 
first one to get the difference 



(A2) 



where ^(e) = 77,+! (e) - 77,(6), and 



,(e) = 1',+i(e)+*P,_i(e)-2T,(e) = ^ [e>F,(6) - 5,,o] . (A3) 



In the last step we have used the original equation of motion 
to simplify the propagator expression. 



If ( A2 1 can be solved for ^A^) for all sites associated with 
the barrier placements, an explicit solution for the probabili- 
ties of all sites m can be obtained. This is the general idea of 
this development. Such a solution can be obtained exactly for 
one or a few barriers through an explicit evaluation of determi- 
nants of manageable size involved in the simultaneous equa- 



tions (A2i. In the case of random and periodic arrangement 
of barrier positions Ps(e), can also be calculated explicitly as 
follows. 



Summation of ( |A2[ l over barrier locations s gives, exactly, 

s s s 

where the new function Jis{e) is defined through 

i"s(e) = Yj ^'■-'■(^'>- 

r 

This /i-function is the sum of propagators among sites placed 
at the left end of the periodically placed barriers. If the barriers 
are placed at random positions, the average can be assumed 
to be independent of the location s, and /ii(e) = iu(e). An 
ensemble average is envisaged here. Then, 



V 1 + ^/"(e) 



is an explicit solution. 

For periodic arrangements of barriers, as used in earlier 
treatments of the problem [18] , the wall location summation 
given by 



(A4) 



can now be done and the result shown to be independent of 
s. As stated in the introduction, the sites to the left of each 
barrier are located at r = Hjl + {H + 1)1 with I varying as an 
integer from -00 to +00. We can thus change the summation 
to s - r = (H + 1)1 and calculate 



(A5) 



To proceed with the evaluation of the sum we use the fact that 
the propagators for chain involve modified Bessel functions. 
It follows that the sum in Eq. ( |A4| i is a geometric sum with 



1 



2F sinh ^ 
tanh (^/2) 



tanh(^(//+ l)/2) 



- 1 



(A6) 



2. Evaluation of Two Sums 

In order to evaluate Eq. (j6]l, we need to perform the follow- 
ing two sums 



(A7) 
(A8) 
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The first can be shown, after some algebra, to be 
sinh(^/2) sinh(^/7) 



/' 1 p) p sinh(^) sinh(^(// + 1 )/2) ' 
whereas the second yields 

r 

H+l sinh(^/2)cosh(^-^g^) 
2Fsinh^ sinh2(«^) 

Now we can use Eqs. \A9) and ( |A10| l in Eq. (|6]) to get 



gH,p{e) p ^ / tanh(f/2) _ , \ 
^ F ltanh(^(H+l)/2) 

(// + 1) sinh(^/2) 



(A9) 



(AlO) 



2F sinh(^) sinh(^(// + l)/2) 
cosh(^(//+l)/2)cosh(^;,) . ' 

sinh(g(//+l)/2) 

H+l-2p I sinh(f /2) sinh(fp) 
2 \F sinh(^) sinh(^(// + l)/2) 



(All) 



Finally, using the expansion of the hyperbolic cotangent leads 
to Eq. dH). 
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